Antiferromagnetic Spinor Condensates are Quantum Rotors 
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We establish a theoretical correspondence between spin-one antiferromagnetic spinor condensates 
in an external magnetic field and quantum rotor models in an external potential. We show that the 
rotor model provides a conceptually clear picture of the possible phases and dynamical regimes of the 
antiferromagnetic condensate. We also show that this mapping simplifies calculations of the conden- 
sate's spectrum and wavefunctions. We use the rotor mapping to describe the different dynamical 
regimes recently observed in 23 Na condensates [l|, Q|. We also suggest a way to experimentally 
observe quantum mechanical effects (collapse and revival) in spinor condensates. 



Bosc-Einstcin condensates occurring in ultracold 
atoms having internal spin degrees of freedom, the so- 
called spinor condensates, offer an exciting addition to 
the family of quantum many body spin systems realiz- 
able in the laboratory Q . Of particular interest are the 
long coherence times and small dissipation rates which al- 
low access to dynamical regimes not available in the solid 
state. Recently there has been considerable experimental 
progress in elucidating the dynamics of spinor conden- 
sates. Such endeavors include dynamics experiments on 
87 Rb atoms for the hyperfme spin-one [!, H| and spin-two 
Q manifolds as well as, most recently, experiments on 
23 Na condensates P, 0- 23 Na spin-one condensates are 
qualitatively different than their 87 Rb counterpart due 
to antiferromagnetic interactions. This leads to ground 
states having zero spin moment as well as disparate dy- 
namical regimes. 

The NIST experiments [l|, Q were performed in a trap- 
ping potential sufficiently tight such that, within a good 
approximation, the bosonic atoms all occupy the same 
spatial mode. This allows the spin dynamics of the sys- 
tem, which is often obscured by spatial variations, to 
be directly probed. The condensate was prepared in an 
initial unstable ferromagnetic state and then allowed to 
evolve freely in time. For small magnetic fields, the sys- 
tem oscillates about the ferromagnetic state, never reach- 
ing zero spin moment (F) = at any time. On the other 
hand, when the magnetic field exceeds a critical value the 
system evolves through (F) = reaching a state point- 
ing in the opposite direction and back periodically, the 
so-called "running phase" trajectories. It was shown that 
these different regimes could be interpreted as being on 
different sides of a separatrix in the phase space of the 
mean-field energy of the system [J Q ■ 

In the single-mode approximation, the full quantum 
Hamiltonian of the system is 
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Here, F = a' a F a paft is the total spin operator where F Q ^ 
are the spin-one matrices and ai,ao,<X-i are bosonic an- 
nihilation operators for each spin state, TV is the total 



particle number, g is the spin-dependent interaction, and 
q is the quadratic Zeeman shift due to an external mag- 
netic field 0- When q = the exact ground state of 
the above Hamiltonian is a condensate of singlet pairs of 
bosons given by Q 
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This ground state is unique and breaks no symmetries. 
However, for large particle numbers, this state becomes 
extremely delicate, being unstable to small external mag- 
netic fields. Thus the observed phases for most experi- 
mental antiferromagnetic systems are more appropriately 
described by symmetry-broken nematic states which are 
well-described by mean-field theory 0, [l(| • This is remi- 
niscent of Anderson's "tower of states" argument for Neel 
ordering in solid state quantum antiferromagnets, despite 
the fact that the true ground state for finite-size bipartite 
lattices can be shown to be a spin singlet [llj . 

In this Letter, we develop a conceptually new approach 
to describe the quantum dynamics of antiferromagnetic 
spinor condensates. In particular, we map the Hamilto- 
nian in Eq. ([T]) onto a quantum rotor Hamiltonian 



Jf = ^L 2 + V(0) 



(3) 



where L is the angular momentum of the rotor, 
/ = Nh 2 /g is the moment of inertia, and V{6) = 

l( N + §) sin2 (0) + sin2 (20) is the external poten- 
tial. The mapping is exact in the sense that the complete 
spectrum of Eq. ([T]) for N bosons precisely agrees with 
the lowest set of eigenvalues of Eq. ([3]) (which has an 
unbounded spectrum from above). A similar procedure 
has been used to derive an exact phase model describing 
bosons in a double- well potential [l2|. One can see that 
the singlet state of paired bosons, Eq. ([2]), corresponds 
to a state where the rotor is delocalized over the entire 
sphere while the symmetry-broken nematic state corre- 
sponds to the rotor being in a position eigenstate. We 
will show how Eq. ([3]) can be used to obtain simple ex- 
pressions for the spectrum and wavefunction of the spinor 
condensate. We then show how the semiclassical limit of 



2 



the rotor system provides a natural interpretation of the 
dynamical regimes of anti-ferromagnetic spinor conden- 
sates observed experimentally Finally, we make 
a prediction to observe quantum mechanical effects (i.e. 
non-mean field effects) in spinor condensates which have 
so far eluded experimental detection. Specifically, we 
show that the abrupt removal of a magnetic field used 
to prepare the system in a nematic state will lead to col- 
lapse and revival dynamics, which cannot be explained 
with mean field theory alone. 

We now proceed with the main technical advance of 
this work: an exact mapping of Eq. (TTJ) onto an effective 
rotor Hamiltonian, thus establishing that antiferromag- 
netic spinor condensates are effective realizations of the 
quantum rotor model. It is most useful to express the 
bosonic creation and annihilation operators as quantities 
that transform as cartesian vectors under rotations. To 
that end we define the operators b x = — (ai — a_i)/\/2, 
by = (ai + a_i)/i\/2, and b z = ao which satisfy bosonic 
commutation relations. It is then straightforward to ex- 
press the Hamiltonian Eq. (fTJ in terms of these operators. 
Specifically, the spin operator is F = — ib^ x b while the 
quadratic Zeeman shift is b\ b z . With these operators we 
construct the complete set of states 
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where SI = (sin(#) cos(<ft), sin(#) sin(0), cos(#)) is a real 
unit vector given by the pair of spherical coordinates 
(9, ft) and N is the number of bosons in the system. For 
simplicity we take N to be even and will comment on the 
odd N case shortly. This wavefunction is the (symmetry 
broken) nematic state pointing along SI. These states 
have the inner product 



(rawing) = (n • si' 



N 



(5) 



Thus, as the number of bosons in the system becomes 
large, states pointing in different directions become or- 
thogonal. 

Interestingly, the spin-singlet state Eq. ((2|) can be con- 
structed by taking equal-weight superpositions of the ne- 
matic state over all directions: 
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as discussed in Refs. 13l. Il4j, This motivates one to use 
the spherical harmonics to construct the orthonormal set 
of states for even £: 
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where ft = 47rN\2 e (^)!/ (^)I(JV + I + 1)! is the 
normalization constant. Such states arc defined for \£\ < 
N, and unless otherwise stated sums for over such states 



are understood to satisfy this restriction. These states 
\£m) can be seen to be eigenstates of the F 2 operator 
with eigenvalue £{£ + !). We finally note that these have 
the following inner product with the nematic states 



(Sl N \£m) = y/f e Y em (n). 



(8) 



With the construction of these two sets of basis states 
|S~2iv) and \£m) in the bosonic Hilbert space we now pro- 
ceed to map the problem onto the rotor Hilbert space. 
This Hilbert space is spanned by the position eigenstates 
| Q) on the unit sphere which arc complete and satisfy the 
orthonormality condition (fi|fi') = <5(S~2 — SI'). These 
states involve angular momentum components for all £ 
and therefore do not suffer the complications that arise 
from Eq. (4) for the |J2jv) states which are only orthog- 
onal in the large N limit . To start we note that a gen- 
eral state in the bosonic Hilbert space can be written as 
a superposition of the spin nematic states with weight 



dCt\n N )ip(n) 



(9) 



We now act with H on this state. If one can find an 
operator H in the rotor Hilbert space such that 
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then a sufficient condition for the time-dependent 
Schrodinger equation (TDSE) in the bosonic Hilbert 
space to be satisfied is the rotor TDSE: = iWt\ip). 
The necessary condition for the rotor model to be a pre- 
cise description for spinor condensates may be less re- 
strictive. 

Our efforts will now be devoted to showing that T-L ex- 
ists and then finding T~L. We consider the two terms of 
the bosonic Hamiltonian Eq. (fTJ) separately. The first 
term, which contains F 2 , is diagonal in the \£m) repre- 
sentation which simplifies the mapping. It is intuitive 
that F a will map to the angular momentum operator in 
the rotor Hilbert space defined 

This can be derived by inserting the completeness rela- 
tions 1 = Y.g m \im){£m\ and 1 = fdn\Sl)(U\ (which act 
in different Hilbert spaces). Using Eq. <j8j) we obtain 

f2 W=E / dnF 2 \£m)^(Ye m \n)(Sl\iP) 
em 

= J dn\n N )(n\L 2 \^) 



where we have used the notation (ri|l& 
Thus we see that 

F 2 -> Kh 2 



(11) 
(12) 
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in the rotor representation. Such a rotor description of 
F 2 was previously noted in 15l4l7|. 

We now move on to mapping the quadratic Zeeman 
term in H to a rotor description. This mapping is more 
complicated since the quadratic Zeeman shift is not diag- 
onal in either the \SIn) or the \£m) representation. Our 
approach will be to express &J& z |f2jv) in terms of |f2jv) 
and its derivatives. Then integration by parts can be 
used to arrive at Eq. (fTU|) . In the analysis we consider 
general quadratic terms of the form b]j)p. We state with- 
out derivation the following identity 



bib p \n N ) = n p (v Q + Jvn a )|n 



N 



(13) 



where V = Ode - 



is the gradient operator on the 



unit sphere. This identity follows from the geometrically 
intuitive relation V Q f2/3 = 8 a /3 — £l a Q/3- We finally note 
that the integration by parts rule for V Q is 



J dn /(n)v„5(«) = [ dsi g{n) [2n a - v a ] /(n). 



this operator commutes with the Hamiltonian Eq. ([3]) 
the constraint imposes no real technical difficulty. The 
case of odd N is similar and is therefore not shown here. 
For this the wavefunction must be antisymmetric and 
the corresponding projection operator running over odd 
£ will also commute with the Hamiltonian. 

We now consider the limiting cases of the rotor Hamil- 
tonian. The simplest situation is when no external mag- 
netic field is present and q = 0. For this the ground 
state is uniformly delocalized over the entire sphere cor- 
responding to the £ = m = spherical harmonic. We 
now consider the case of small magnetic field such that 
g ^> q > 0. For this case the first term in the rotor 
potential V{9) dominates and serves to localize the ro- 
tor about the poles. In this limit, we can expand the 
potential to quadratic order about the 9 = minimum 
and the Hamiltonian becomes that of a two-dimensional 
harmonic oscillator [l8|. The spectrum for the lowest 
energies are then 



(14) 



(19) 



Using Eqns. ©, ([I3"j). and (T3]) we obtain 
&LM*} = j dQip(£l)Q/3 (V a + NCl a ) \CIn) (15) 



(for even n with multiplicity 2n+l) and the ground state 
wavefunction is 
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2qN 2 
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That the 



From this we can read off the equivalent operator acting 
in the rotor space which corresponds to b^bp: 



blbp -> (N + 3)Q a flp - ^/?V Q - 6 a p. 



(16) 



Using the mappings in (fT2j) and (fT6|) restricted to the 
case a = /3 = z, we finally arrive at the operator %: 
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where V z = — sm(9)dg and we have dropped a constant 
term. While W has a real spectrum, it is not Hcrmi- 
tian. It is therefore advantageous to apply a similarity 
transformation to render it Hcrmitian. Defining 



,ye = e F He- F 



(18) 



with F = — ^ cos(2#) we arrive at Eq. and the map- 
ping is complete. We note that with this transformation, 
the wavefunctions ip(£l) governed by Ji?, when entering 
Eq. ([9]) must be accompanied by a factor of e~ F . 

This equation is the model for a quantum rotor under 
an external potential. Since we are taking the case of even 
N the wavefunctions must satisfy the constraint ip(fl) = 
This condition can be interpreted as constrain- 
ing the ends of the rotor to be bosonic particles, requiring 
the rotor wavefunction to be symmetrical under their in- 
terchange. This constraint can be enforced with the pro- 
jection operator V = ^ even ^ Y^ m =-l \ Y tm)(Y lm \. Since 



energy states are evenly spaced and have the spectrum 
given by Eq. (|19p in this regime is not immediately clear 
from a direct analysis of the original bosonic Hamiltonian 
Eq. ([T]). In order for this harmonic oscillator description 
to be valid we must have the condition 9 <C 1. Away 
from this limit the rotor will delocalize and approach the 
singlet state. For a large particle number N we therefore 
see that any small external magnetic field will tend to 
drive the system to the symmetry broken nematic state 
as described by the mean field theory @,[l3l- For higher 
magnetic field we see that when q > 2g a local minimum 
appears along the equator 9 = | though the global mini- 
mum will remain at 8 = 0. This leads to stationary states 
localized about the equator. Such states are analogous 
to the "7r-states" occurring for a scalar condensate in a 
double- well potential [l9j]. However, as in the double- well 
case, transforming this wavefunction back to the bosonic 
Hilbcrt space can significantly alter its structure 12 1 . 

Having described quantum mechanical states of Eq. 
in various limiting cases we now proceed to a semi- 
classical analysis of its dynamics which is relevant to the 
recent experimental results [H, 0- The Lagrangian de- 
scribing the motion in the semiclassical limit is 



C = -I [e 2 + sin 2 (6>)</> 2 ) - V{9). 
The equation of motion for this is 

fe = i--fl P 2 - 

sin 3 (0) 
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FIG. 1: Collapse and revival of (ip(t)\Q, 2 z \ip(t)} starting from 
a state localized about poles with width 9 = 0.1. 



where 



sin 2 (9)<j) is a constant of motion. As before 



we start by considering the limiting case g ^ q > 0. For 
this case we can drop the second term in V . The first type 
of motion we consider is when the rotor remains close to 
the minimum at the poles at all times. The potential can 
then be expanded to quadratic order in 9 and analytic 
solutions can be found. One solution is where the rotor 
oscillates through the poles: 9(t) — 9ocos(ujt), = 0. 
Another solution is where the rotor precesses about the 
poles: 9 = 0, cj)(t) = ut. Both of these solutions have the 
cigcnfrequcncy to = y/2gq/Ti which corresponds to the en- 
ergy scale appearing in the spectrum from the quantum 
mechanical analysis Eq. ([11?]) . The second type of motion 
we consider is where the rotor has enough energy to over- 
come the potential barrier near the equator and explore 
both hemispheres in its trajectory. These are precisely 
the oscillating phase solutions experimentally observed 
in [H, Q ■ Finally, a third type of motion is possible when 
q > 2g. As described above, for this case there is a lo- 
cal minimum at the equator. Therefore for this situation 
there will be trajectories which remain localized about 
the equator. 

We now apply the rotor description to the quantum 
dynamics of antifcrromagnetic condensates in the single 
mode regime, which is known to manifest rich behavior 
[20l - l22| . Here we consider preparing the system in the 
symmetry- broken nematic state given by Eq. (|20[) . and 
then rapidly turning the magnetic field off and allow- 
ing the state to evolve freely. We note that according to 
the semiclassical theory (or by using the Gross-Pitacvskii 
equation) the nematic wavefunction will remain at the 
pole and not evolve temporally. The quantum mechan- 
ical dynamics, however, is markedly different. By dy- 
namically evolving the wavefunction, Eq. (|20[) . with the 
quantum rotor Hamiltonian Eq. ([3]) with q = 0, it can 
be seen that the state will undergo periodic collapse and 
revival at the characteristic frequency h/I. For instance, 
provided the initial state is sufficiently localized 9 <C 1, 
one can show that 



The evolution of this function over a single period is 
plotted in Fig. [T] The localized nematic state rapidly 
collapses to states with substantial weight contributions 
from other regions of the unit sphere, and then fully re- 
vives at the end of the period. By applying the Poisson 
resummation formula to Eq. ([25)1 it can be seen that the 
evolution is a train of localized pulses separated by a 
fourth of the time period. This behavior can be directly 
seen experimentally by measuring the time dependence 
of (ajao) after the turning off the magnetic used to pre- 
pare the system in the polar state. We note that since 
the magnetic field couples only to the spin degrees of free- 
dom, the above procedure will not excite spatial modes 
of the condensate for sufficiently tight traps. The quan- 
tum collapse and revival of Fig.[TJis a direct consequence 
of the rotor mapping of spinor condensates. 

In conclusion, we have established a correspondence 
between antiferromagnetic spinor condensates and quan- 
tum rotors. We have shown that this mapping offers a 
considerable conceptual as well as technical advance in 
understanding the properties of spinor condensates. We 
use the mapping to address recent experimental results 
0, [H and to analytically predict a collapse and revival 
process (which is a direct experimental signature of quan- 
tum effects). We point out that it should be possible to 
provide similar quantum rotor mappings for condensates 
with larger spin. 
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